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Using a fully self-consistent quantum statistical model, we demonstrate the possibility of Bose condensation 
of nuclei in heavy ion collisions. The most favorable conditions of high densities and low temperatures are 
usually associated with astrophysical processes and may be difficult to achieve in heavy ion collisions. None- 
theless, some suggestions for the possible experimental verification of the existence of this phenomenon are 
made. 
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For some time, there has been an intrinsic interest in Bose 
condensation of particles with a concomitant phase transition 
in nuclear physics. Studies of multifragment final states over 
a broad range of incident energies may provide information 
about liquid-vapor and other types of nuclear phase transi- 
tions, and ultimately information about the equation of state 
of nuclear matter [1-12], provided that multifragment disin- 
tegrations are related to bulk instabilities of nuclear matter at 
low densities. Most experimental observations [8] provide 
end-point information and do not necessarily provide infor- 
mation about transient phenomena. Here we demonstrate the 
possibility that Bose condensation of nuclei is such a tran- 
sient phenomena, which has thus far been unobserved. The 
aim of this Rapid Communication is twofold: first, to estab- 
lish the possibility of Bose condensation of nuclei and the 
concomitant phase transition in nuclear physics; second, to 
suggest some specific experimental observations for its veri- 
fication. The smallness of the finite volumes being consid- 
ered could mask a clear signal of the phase transition. How- 
ever, by exercising sufficient care, experimental signatures 
for the onset of Bose condensation may be observable. 

Bose condensation of pions was predicted earlier [6] to 
occur in the expansion phases of nuclear collisions. There is, 
however, a subtle difference between Bose condensation of 
nuclei and that of pions. In the formation of nuclei, one 
requires several nucleons to be present in close physical and 
momentum proximity. Therefore, matter present must have 
some minimum density for condensation to occur. On the 
other hand, pion production is basically a one-nucleon phe- 
nomenon and therefore does not depend upon the presence of 
a minimum density of nuclear matter. Hence, pion conden- 
sation is observable in end-point measurements, such as total 
yields or cross sections, but condensation of nuclei is not. Its 
observation must be related to pre-ultimate stages of rarefac- 
tion (e.g., compression). 

To focus on the possible existence of the phenomenon, a 
simple model is assumed, which ignores interactions be- 
tween the particles. Such a model should be reasonable for 
the low densities and temperatures considered herein. More 
detailed and reliable predictions would certainly require a 
more sophisticated and realistic model; nevertheless, a 
simple model will suffice to demonstrate the possibility of 
observing this phenomenon. 
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We begin by considering an initial assembly of N {) neu- 
trons and Z () protons which have been compressed and are 
evolving in thermodynamic equilibrium at an excitation en- 
ergy (temperature) T and density p to form light ions: 12 C, 
"B, 10 B, 9 Be, 7 Li, fi Li, 4 He(a), 3 He(/t), ? H (/), 2 H (d), 
together with neutrons (n) and protons ( p ). (Although the 
existence of a fully equilibrated system in heavy ion colli- 
sions is still debated, it is generally believed [4,5,9-12] that 
equilibrium thermodynamics is a viable theory which gives 
reliable predictions.) In the model, charge and baryon num- 
ber are conserved by requiring that 

2 ZiNj = Z„, 


2 HiNi = N 0 , (I) 

i 

where the N t refer to the number of particles of species / .In 
equilibrium, the chemical potential is 

Mr = z ,7*p + «,/*« + * : M (2) 

where z i and n i are the number of protons and neutrons in 
the /th species, and e, is its binding energy. These particles 
have finite dimensions and initially move within a volume 
V {) . For convenience, we treat the particles as point particles 
(no physical dimensions) which move in a reduced volume 
given by 

v^v,, -Saw, O) 

i 

where V i is the eigenvolume of the /th particle. The number 
(matter) density is 

P = (W 0 + Z 0 )/VV (4) 

The density of point particles is 

P pt =(yv 0 +z 0 )/v pt , (5) 
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where p p{ >p since V pt < V () . Our results will be presented in 
terms of p pt . 

For a distribution of fermions, quantum statistics [13] 
yields 

Ni — 2gjV p{ Tr- l/2 X~ 3 F f D (fij/k B T), (6) 

where g i = 2S i + 1 is the spin degeneracy, ju f is the chemical 

potential, is the thermal wavelength for the ith particle of 
mass m l , 

\ ( = 2Trfi(2'nm i k B T) 1/2 , (7) 

where k B is Boltzmann’s constant, and the species index is i 
(=/?, /z, /z, /, 7 Li, 9 Be, 11 B) as appropriate. Values for the 
Fermi-Dirac integral functions 

r x ,/2 

Ffd(>)= J 0 dx > (8) 

are tabulated in the literature [14] and are increasing func- 
tions of v. For bosons [13,15,16], we have 

AE= [exp( - fJL,lk B T) - 1 ]" 1 +g ( Vpt^i V BE ( - //,,*/ k B T), 

(9) 

where the first term on the right-hand side of Eq. (9) gives 
the number of condensed bosons, the second term is the 
number of noncondensed bosons, and the species index is i 
( = d, a, 6 Li, 10 B, I2 C). In Eq. (9), the Bose-Einstein inte- 
gral functions [15] are given by 

f* * l/2 

Fbe(v)=\ zrnr-rdx. ( 10 ) 

Jo e “i 

The evolution of the system is studied by self-consistently 
solving Eqs. (6) and (9) while simultaneously requiring Eqs. 
(1) and (2) to be satisfied, until a convergent solution has 
been found for a specified temperature and density p pi . The 
density is then increased for that temperature until another 
convergent, self-consistent solution is found. Using this pro- 
cedure, we find that there is a maximum density, for each 
temperature, beyond which no self-consistent solution exists. 
This occurs when the magnitude of the chemical potential for 
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FIG. 2. Percentage of condensed bosons as a function of density 
ratio at temperatures of 2.5 and 5 MeV. 

one of the bosonic species (usually a or 12 C) becomes small 
leading to Bose condensation and a concomitant phase tran- 
sition. 

To illustrate the application of these methods, consider a 
model system, of A o = 40 neutrons and Z 0 = 40 protons, 
such as might exist in a 4() Ca- 4n Ca central collision. Figure 1 
displays particle distributions for 7=2.5 and 5 MeV as a 
function of the density ratio p pl /po where p n = 0.17 fm"\ 
Note that the ion distributions vary substantially with 7 or p. 
From these results it is clear that bosons are more abundant 
than fermions. The most abundant ion species are l2 C, 
*°B, and 4 He. Figure 2 displays the percentage of condensed 
particles for 7=2.5 and 5 MeV From these results we note 
that there is a finite probability of producing condensed alpha 
particles and condensed 12 C nuclei at densities occurring in 
heavy ion collisions (p^2p 0 ), and even for densities below 
that of normal nuclear matter (p<p n ). Also, the percentage 
of condensed bosons increases with increasing density. Fig- 
ure 3 displays the chemical potentials for the most abundant 
bosons. All other boson species have values which are 
smaller (more negative) than those displayed in the figure. 

From the results of these analyses with a simple, yet re- 
alistic physical model, it is apparent that Bose condensation 
of nuclei, such as 4 He or 12 C, is a real possibility in central 
collisions of heavy ions. 

The question is, how do we observe it? Although the most 
favorable conditions (for producing condensed bosons) of 
high densities and low temperatures (preferably less than 8 
MeV) may be difficult to achieve in laboratory experiments, 
the reasonably large percentage of condensed particles, indi- 
cated by our calculations, may provide a measurable signa- 
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FIG. 1. Distribution of most abundant particles up to densities of 
2.5p () . Curves for temperatures of 2.5 and 5 MeV are displayed. 
Other species have abundances which are smaller than the lowest 
curve shown here. 


FIG. 3. Chemical potential as a function of density ratio for the 
most abundant bosons. Curves for temperatures of 2.5 and 5 MeV 
are displayed. Other bosons have values smaller than those dis- 
played here. 
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ture in more realistic laboratory conditions. The most com- 
mon observable for any distribution of particles is to measure 
total yields or cross sections. A measurement of the total 
yields, however, would not confirm or deny the existence of 
condensed bosons. From Fig. 1, we note that the total yield 
of bosons (condensed and noncondensed) varies with tem- 
perature and density. From Fig. 2, the relative proportions of 
condensed to noncondensed bosons also varies. Conse- 
quently, experimental observation or detection must be based 
upon methods which distinguish between these two catego- 
ries of bosons. One such method of detection might involve 
measuring the momentum distributions of a boson species 
(e.g., 12 C or 1 2 3 4 5 6 7 8 He). In the center-of-mass system, the con- 
densed bosons will be produced with near-zero momenta, 
whereas the noncondensed bosons will have much larger mo- 
mentum values. Therefore, finding bosons with near-zero 
momenta in the center-of-mass system would be a clear sig- 
nal that Bose condensation has occurred. In the laboratory 
frame, these condensed bosons would be detected at 0° to 
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the beam with a momentum equal to that of the collision 
system center of mass. There should be no significant diffi- 
culty in detecting them [17]. 

Finally, we emphasize that the findings and conclusions of 
the present work, although definitive in nature, deserve fur- 
ther study and refinement using more sophisticated models 
than the one used here. Nevertheless, the specific signature 
indicates that experimental verification of our proposal 
should be possible even though difficult because of the finite 
sizes of nuclear ensembles. In view of the importance of the 
phenomena to studies in nuclear physics and astrophysics, 
verification of its existence or nonexistence merits concerted 
efforts by experimental investigators. 
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